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State University of New York at Albany, Albany, New York 12222 
(Dated: February 7, 2008) 

Abstract 

Using a 13.7 fb^"*^ data sample collected with the CLEO II and II. V detectors, we report new 
branching fraction measurements for two Cabibbo-suppressed decay modes of the meson: 
B{D+ 7r+7r°) = (1.31 ± 0.17 ± 0.09 ± 0.09) x 10"=^ and B{D+ K+K^) = (5.24 ± 0.43 ± 
0.20 lb 0.34) X 10~^ which are significant improvements over past measurements. The errors reflect 
statistical and systematical uncertainties as well as the uncertainty in the absolute branching 
fraction scale. We also set the first 90% confidence level upper limit on the branching fraction of 
the doubly Cabibbo-suppressed decay mode B{D~^ — >■ K'^tt^) < 4.2 x 10"*^. 
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To lowest order, weak decays of mesons may be described by the six quark- diagrams 
shown in Fig. Q external W-emission, internal W-emission, W-exchange, W-annihilation, 
horizontal W-loop, and vertical W-loop 0. When using these diagrams to describe pro- 
cesses, dynamical assumptions are often made regarding the relative size of their amplitudes 
as well as the nature of interference terms between diagrams. Measurements of hadronic 
decays of mesons give insights into these assumptions as well as new information on 
the violation of SU{?>) flavor symmetry {SU{?>)f) violation, isospin symmetry, and doubly 
Cabibbo-suppressed decays. 

SU{'i)F symmetry breaking is of current interest because of — mixing studies; it 
has been shown that the mass and width differences (x, y) of the CP-eigenstates of neutral D 
mesons can be generated by the second order 5'f/(3)i? symmetry breaking 2]. Understanding 
the size of these effects may be important to unravel any non-Standard Model contributions 
to D^D^ mixing. Such understanding is only possible if SU (3) p violating effects are well- 
determined. We report new measurements of the decay modes tt+tt^ and — »• 
K^gK^ , which are useful for the estimation of SU{?))f violating effects in the D meson 
system. 

Predictions based on isospin symmetry are generally considered to be more reliable than 
SU{?>)f predictions because of the near degeneracy in mass of the u and d quarks. Using 
measurements from this analysis as well as data from the Particle Date Group (PDG) j^, 
we determine the isospin amplitudes and phases for the D — > tttt system. 

Doubly Cabibbo-suppressed decays (DCSD) of charm mesons involve c ^ d and s ^ u 
quark transitions. Currently, there are only four measured DCSD decay modes j^. Measure- 
ments of such modes will lead to improved understanding of SU (3) f and other Standard 
Model predictions. Such modes are also important for neutral D-mixing measurements, 
where a significant background is from DCSD decays. In this paper we report the first 
upper limit on the branching fraction of the DCSD decay K'^tc^. 

This analysis uses data collected with two configurations of the CLEO detector at the 
Cornell Electron Storage Ring (CESR): CLEO II || and CLEO II. V Q. The total inte- 
grated luminosity of the data sample is 13.7 fb~^. The CLEO detector is a general purpose 
spectrometer with excellent charged particle and electromagnetic shower energy detection. 
In CLEO II the momenta of charged particles are measured with three concentric drift cham- 
bers between 5 and 90 cm from the e+e~ interaction point. In the CLEO II. V configuration 
the innermost drift chamber was replaced by a 3 layer silicon vertex detector. Charged 
particles are identified by means of specific ionization measurements (dE/dx) in the main 
drift chamber. The tracking system is surrounded by a scintillation time-of-fiight system 
and a CsI(Tl) electromagnetic calorimeter. These detectors are located inside a 1.5 T su- 
perconducting solenoid, surrounded by an iron return yoke instrumented with proportional 
tube chambers for muon identification. 

Charged pion and kaon candidates were required to pass minimum track-quality criteria. 
Kaon (pion) candidates had to have a specific ionization within two (three) standard devi- 
ations (a) of that expected for a true kaon (pion). We combined pairs of electromagnetic 
showers in the calorimeter to create 7r° candidates. Candidates with a reconstructed mass 
within 2.5 a of the nominal vr" mass were kept for further studies. We obtain Kg candidates 
by reconstructing the decay mode Kg —>■ tt+tt^. We required daughter tracks to have an 
impact parameter in the plane transverse to the beam greater than three times the measure- 
ment uncertainty and that the probability of the returned from the vertex fit for pairs 
of daughter tracks was required to be greater than 0.001. Kg candidates also had to have a 
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FIG. 1: Six lowest order quark diagrams for a meson decaying into two mesons f])|: (a) external 
W-emission, (b) internal W-emission, (c) W-exchange, (d) W-annihilation, (e) horizontal W-loop, 
(f) vertical W-loop. Dashed lines represent W boson. 



reconstructed mass within 3.0 a of the nominal mass. 

In order to reduce backgrounds, we required that candidates come from the decay 
D*^ D^TT^, with the mass difference (AM) of the reconstructed D*^ and to be 
within 2.5 a of the known value 0. We required all D*^ candidates to have a normalized 



momentum {xd* = Ipd* l/y (5/2)2 — mj^^^ greater than 0.6 and all candidates to have 
a cosOfi value between ±0.8, where 9h is the angle between the charged daughter track in 
the rest frame of the and the reconstructed momentum vector in the rest frame of 
the D*~^ meson. To insure that we obtained only one candidate per event, we selected 
candidates with the lowest value for 

, _ (AM - AMpDG? ~ "^yy ^^ 

where i indexes the 7r° candidates in this decay. Given the large uncertainties in absolute 
branching fractions we present our results as ratios of the branching fraction of the decay 
mode under study to that of a normalization mode: — A'~7r"'"7r+ for n~^7r^, K'^ir'^ 

and D+ -> K'^s7i+ for D+ K^s^^ . 

To extract the yield for each mode, we performed an unbinned maximum likelihood fit 
for two components (signal and background) using the following observables: m/), the mass 
of the reconstructed meson, the normalized momentum of the D*^ meson, Xd*-, and 
cos 9h- Using a GEANT-based simulation of the CLEO detector as well as sideband data 
we determined probability density functions (PDF) for each observable describing the shape 
of the data for signal and background events for each decay mode. The probability that a 
candidate is consistent with signal or background is given by the product of these PDFs. The 
likelihood is given as the product of these probabilities over all candidates; maximization 
of the log of the likelihood gives us the signal and background yields. Projections of the 
likelihood fit to the mass for our three decay modes are shown in Fig. |21 Using simulated 
signal and background events, we measure the efficiency of our analysis method for each 
mode, enabling us to determine the total number of signal events in our data sample for 
each decay mode. Table HI lists raw yields and efficiencies for all decay modes. 



4 



> 

CD 



60 



1 1 60603-002 



40 



c20 

CD 

m 



1 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 
( Pi) 


1 1 1 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 



> 

CD 



1.77 
60 



CD 40 



I 20 

CD 

iS 



1.82 1.87 1.92 

M(Jt+jrO) (GeV) 

1 — \ — \ — \ — \ — \ — \ — \ — I — \ — \ — \ — I — \ — r 

(b) 



1.97 




* I I I I I I I I I I* I ^v^*^ 



> 





1.83 1.85 1.87 1.89 1.91 
M(K+K°) (GeV) 

80, 



CD 

£240 
(/) 

CD 

m 



III 


1 1 1 1 1 1 1 1 1 

(c) - 


1 1 1 1 


1 1 1 1 1 1 1 1 1 1 1 1 



1.75 1.79 1.83 1.87 1.91 
M(K+jt°) (GeV) 

FIG. 2: Invariant mass distributions for (a) D+ 7r+7r°, (b) D+ K+K'^, and (c) D+ K+tt^ 
candidates. The points represent the data and the hnes are the projections from the maximum 
hkehhood fit. 



We considered systematic uncertainties from experimental resolution, efficiency determi- 
nation, and PDF parameterization. The first two contributions are small and the systematic 
errors are dominated by uncertainties in the PDF parameterization. We studied this sys- 
tematic effect for each mode by simultaneously modifying every PDF parameter within its 
uncertainty. We extracted the yield from the data after each modification to produce a 
distribution of yields. We defined the systematic uncertainty due to PDF parameterization 
as the 68% limits for these distributions. 

Combining the systematic error study with the yields and efficiencies given in Table Ewe 
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TABLE I: Yields from the maximum likelihood fit with statistical errors and reconstruction effi- 
ciencies. 



Mode 


Yield 


Efficiency 




171.3 ±22.1 


(6.20 ±0.11)% 




277.7 ±20.8 


(4.94 ± 0.23)% 




34.3 ± 20.9 


(6.08 ± 0.22)% 




12898.0 ± 156.6 (6.74 ± 0.12)% 




1434.7 ±48.0 


(4.83 ± 0.23)% 



obtain the following results: 

}3{D+ tt+ttO) 



B(D+ ^ K-7r+7i+) 

B{D+ vr+ETO) 
B{D+ K+Tc^) 



0.0144 ±0.0019 ±0.0010 
0.1892 ±0.0155 ±0.0073 
0.0029 ±0.0018 ±0.0009 



where the first error is statistical and the second error is systematic. The results supersede 
previous CLEO measurements 0, 1^ . 

In order to determine the absolute branching fractions, we combine our results with the 
PDG values oiB{D+ K--K+-K+) = (9.1 ±0.6)% and B{D+ n+K^) = (2.77±0.18)% 
and find 

i3(L>+ ^ 7r+7r°) = (1.31 ± 0.17 ± 0.09 ± 0.09) X 10~^ 
B{D^^K^K^) = (5.24 ± 0.43 ± 0.20 ± 0.34) X 10"^ 
B{D^ ^ K+TT^) = (2.64 ± 1.64 ±0.82 ±0.17) X 10"^ 

where the third listed uncertainty comes from the error of branching fractions of the nor- 
malization modes. 

With no significant signal being observed for the doubly Cabibbo-suppressed decay 
—>■ K^Tx^ we determined the 90% confidence level upper limit for this branching frac- 
tion. Our method for obtaining the upper limit involved creating 1000 new data sets with 
the same number of signal and background events as our data sample. In order to include 
systematic uncertainties in our upper limit, we also modified the PDF parameters in the 
manner described for our branching fraction calculation. Using this method, our upper limit 
is 

B{D+ K+n^) < 4.2 X 10'^ at 90% C.L. 
In the limit of SU{3)f, the following ratio is expected to be unity 



cd 



where the Vcs and Vcd arise because of the different quark transitions in the two decays and 



the factor of 2 arises because of the \/l/2 term in the normalization of the vr'' wavefunction. 
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Using = 4.45 ± 0.32 j^, the yields and efficiencies (Table HJ obtained from our 

analysis and combining statistical and systematical uncertainties in quadrature, we find 

Ri = 1.84 ±0.38 

slightly inconsistent with theoretical expectations that SU{3)p symmetry breaking effects 
are about 30%. 

It is believed that in the D meson system the interference between external and internal 
W-emission decay amplitudes (Figure H} is destructive. In order to test this assumption we 
calculate the ratio 

1 T{D+ K+K"^) _ T{D+ K+K'^s) 
^ ~ 2 ^ T{D+ tt+ttO) ~ T{D+ n+n^) ' 

which in case of destructive interference should be greater than 1. Besides a small contri- 
bution from the W-annihilation diagram the decay in the numerator, K^Kg, can 
be described using an external W-emission diagram. Whereas both the external and the 
internal W-emission amplitudes contribute to the decay in the denominator, — >■ tt+tt^. 
Experimentally, we find using our yields and efficiencies from Table |l] 

i?2 = 2.03 ± 0.32 

indicating that the interference between external and internal W-emission is indeed destruc- 
tive. 

Final state interactions (FSI) are significant in charm decays. Using our measurement 
for — > TT+TT^ and the PDG values for D'^ it~^tt~ , tt^tt^ Q we can gain some insights on 
these effects by determining isospin amplitudes and phases for the D ^ tttt system. The vrvr 
final state may have an isospin value of or 2. Writing the amplitudes for the 1 = state 
as Aq and the 1 = 2 state as A2, we obtain the following relation: 

2 p+O 

= |(r+- + foo) - r+o 

where F*^* = r{D^ vr'^vr^) and a, b represent the charges of the pions. Since isospin 
amplitudes are complex, measuring the phase between them is necessary to obtain full 
information about the amplitudes. The phase is written as 

3F+- - 6F°o + 2F+° 
coso = — _ 

4y2r+o^|(F+- + Foo) - F+o 

We find IA2M0I = 0.421 ± 0.044 and cos 5 = 0.042 ± 0.195. These results supersede a 
previous CLEO measurement |7|. The large relative phase between the isospin amplitudes 
indicate that there are significant FSI effects in the D — >■ vrvr system confirming our earlier 
results A similar observation has been made recently by the FOCUS collaboration [,10| . 

In summary, we have obtained measurements for two singly Cabibbo-suppressed 
decay modes: 8{D+ tt+tc^) = (1.31 ± 0.17 ± 0.09 ± 0.09) x 10-^ and B{D+ K+K^) = 
(5.24 ± 0.43 ± 0.20 ± 0.34) x 10~^. We also present an upper limit on the DCSD mode 
B{D^ —>■ K'^n^) < 4.2 x 10~^ at the 90% CL. Our experimental measurements confirm the 
destructive nature of the interference term between the external and internal W-emission 
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diagrams and indicate significant SU{3)f symmetry breaking. An isospin analysis shows 
that FSI effects are important for hadronic decays of D mesons. 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent 
luminosity and running conditions. This work was supported by the National Science Foun- 
dation, the U.S. Department of Energy, the Research Corporation, and the Texas Advanced 
Research Program. 
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